A Limit on Primordial Small-Scale Magnetic Fields from CMB Distortions 
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Spatially varying primordial magnetic fields may be efficiently dissipated prior to the epoch of 
recombination due to the large viscosity of the baryon-photon fluid. We show that this dissipation 
may result in observable chemical potential jj, and Compton y distortions in the cosmic microwave 
background (CMB) spectrum. Current upper limits on fi and y from FIRAS constrain magnetic 
fields to have strength Bo < 3 x 10^* Gauss (scaled to the present) between comoving coherence 
length ft; 400 pc and ft; 0.6 Mpc. These represent the strongest upper limits on small-scale primordial 
magnetic fields to date. 
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At present not much is known about the existence or 
absence of magnetic fields in the early universe. There 
is a plethora of proposed scenarios for the generation of 
primordial mag netic fields Most of these scenarios 

operate by "causal" mechanisms, such that the coherence 
length of the resulting magnetic field is usually limited by 
the horizon scale at the epoch of magnetogenesis (e.g., ~ 
1 pc at the QCD transition), unless magnetic field energy 
density is subsequently transferred from smaller to larger 
scales by inverse cascades ||]. 

Limits on the total magnetic field energy density ex- 
tant at the epoch of Big Bang nucleosynthesis (BBN) 
may be derived (i.e. Bq ^ 10^^ Gauss) from the pos- 
sible increase of ^He production beyond observational 
bounds (Here, and henceforth, we quote magnetic 
field strengths at the values they would have at the 
present epoch, i.e. after cosmological expansion.) Never- 
theless, such strong magnetic fields usually do not survive 
neutrino induced dissipation prior to BBN rendering 
this limit useless in most cases. Fields on large scales may 
be limited by the cosmic microwave background (CMB) 
and by Faraday rotation measures of light from high red- 
shift quasars. A homogeneous cosmic magnetic field on 
the scale of the present horizon is limited by its effect on 
the CMB due to an induced anisotropic cosmic expansion 
to fall below Bq ^ 3x 10~^ Gauss [D while the Faraday ro- 
tation limits are slightly weaker at i?o 10~^ Gauss [1). 
For smaller coherence lengths, Faraday rotation mea- 
sures constrain primordial fields with a 50 Mpc coherence 
length at present to have i?50Mpc ^ 6 x 10~^ Gauss and 
-Bmpc ^ 10^* Gauss for Mpc coherence lengths [||. In 
addition, a spectrum dependent limit on stochastic mag- 
netic fields of Bq ^ 10"^ Gauss on ^ 0.1 Mpc scales has 
been quoted due to the possible CMB distortions from 
magnetic field generated gravitational waves llfl] . 



Future high-precision CMB satellite observations 
should be able to detect (or limit) large-scale mag- 
netic fields. Observational proposals include the de- 
tection of rotational (vector) perturbations below |Tl| ] 
and above the Silk damping scale, a depolarization 
of the cosmic microwave background by Faraday rota- 
tion as well as distortions of the acoustic peaks 
due to a changed (magneto-) sonic propagation veloc- 
ity at the epoch of recombination These effects 
may be observable for magnetic field strengths above 
Bo ^ 3 X 10~^ — 10~^ Gauss given anticipated sensitivities 
of the satellite missions, and under the assumption that 
magnetic field coherence lengths are not much smaller 
than the comoving Silk damping scale (~ lOMpc). 

In this letter, we impose a new and fairly generic limit 
on pre-recombination era magnetic fields on small scales 
based on stringent limits on the deviation of the CMB 
spectrum from a perfect blackbody. 

We have recently shown that spatially varying, 
stochastic, magnetic fields undergo significant dissipa- 
tion before recombination beyond that of magnetic dif- 
fusion 01 . To see how fields dissipate, imagine initial 
conditions of a randomly magnetized baryon-photon fluid 
with zero peculiar flows and at uniform fluid density. The 
non-force-free component of the field accelerates the fluid 
resulting in the conversion of magnetic energy into ki- 
netic energy of the fiuid. For all length scales. A, with 
magnetic relaxation time, Tr X/va shorter than the 
Hubble time, Tr < (where H is the Hubble pa- 

rameter), the resulting fluid motions may be described 
by a superposition of Alfven- and slow magnctosonic- 
waves with approximate energy equipartition between 
magnetic and kinetic energy. Here va is the Alfven ve- 
locity, which well before recombination (and on scales 
larger than the photon mean free path) is given by 
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VA = B/^Att{p^+p^) = 3.8 X IQ-'^c (Bq/IO-^G), where 
/9-y and p-y are photon- energy density and pressure, re- 
spectively. (In what follows we adopt natural units with 
the speed of light c~l.) The generated kinetic energy is 
efficiently dissipated into heat by the considerable shear 
viscosity of the pre-recombination fluid, 77, which arises 
from the long photon mean free path. The dissipation of 
energy drawn from the initial magnetic field resumes until 
the magnetic field reaches a force- free state, i.e., does not 
exert any magnetic stresses on the fluid. For most mod- 
els of magnetic field generation in the early universe, if 
an initially random magnetic field relaxes to a final force- 
free field configuration (for example, straight field lines), 
the final state will only have a negligible fraction of the 
initial magnetic field energy. 

A detailed calculation of the dissipation times for 
Alfven, slow- and fast- magnetosonic perturbations on a 
uniform magnetic field background as a function of mode 
wavelength A was performed in Q . Similar to the damp- 
ing of baryon-photon sound waves by photon- shear vis- 
cosity and heat conductivity (i.e.. Silk damping), effective 
damping of magneto-hydrodynamic perturbations occurs 
in the photon diffusion limit, when the photon mean free 
path <C A. At a given cosmic time, t, modes with A 
smaller than the photon diffusion length, ~ (^^i)^/^, 
are efficiently dissipated. It should be noted that al- 
though the dissipation scales have been computed in |^] 
in the linear regime, it seems difficult to avoid dissipation 
by non-linear effects unless power on small scales can be 
efficiently transferred to large scales increasing the mag- 
netic field coherence length faster than the increase in 
the dissipation scale. Note also that a subset of the full 
non-linear MHD equations exhibits the same damping as 
linear modes p5[ |. 

In general, damping of Alfven and slow magnetosonic 
waves is significantly different from damping of sound and 
fast magnetosonic waves when 1;^^ <C 1, corresponding to 
present magnetic field strengths Bo <C 3 x 10^® Gauss. 
When the magnetic field relaxation time, ^ X/va, is 
longer than the dissipation time scale (by photon diffu- 
sion), Td ~ A^/Z-y, the magnetic field configuration be- 
comes overdamped. In this limit, modes take longer 
than Tr to accelerate the fluid since the fluid is strongly 
"dragged" . For some modes, this overdamping occurs 
while still in the photon diffusion limit {1^ <C A). Such 
modes do not dissipate until later when the photons are 
free-streaming on scales, l-y ^ A. In this case, the ef- 
fective damping scale is much shorter than the photon 
diffusion length, ^ VAd-y (cf. 0). 

In this letter, we are mostly concerned with damp- 
ing of magnetic flelds in the photon diffusion limit, since 
this leads to the strongest constraints on small scales. 
In the dissipation times were computed in the WKB 
approximation, i.e., with mode amplitudes evolving as 
bk{t) — &fc(io)exp(i J^^Ludt), where to is mode frequency 



and k is comoving wavevector. In this expression b de- 
notes a "comoving" magnetic fleld strength, scaled to 
the present epoch, with b = ba^, where b is the red- 
shift dependent magnetic fleld and a is the scale factor 
(a = 1 at the present epoch). Dissipation of magnetic 
flelds manifests itself through an imaginary part in the 
mode frequency. The dispersion relations for Alfven and 
slow magnetosonic waves in the photon diffusion limit are 
identical to lowest order in b 

In this expression 3(p^ -I- p^)??' — t], where 77 is shear 
viscosity, p-^ and as above, and R = 2>pi,/{^p^) « at 
redshifts z ^ 10"*, with pb baryon density. Photon shear 
viscosity is given by 77 = (4/15)^^/7 with the photon 
mean free path l~f = 1/ {ne<JTh) determined by Thomson 
scattering on electrons. 

One may Fourier transform the stochastic mag- 
netic ^fleld by writing, b{x) = {V^'^ / {2TTkcf''^) 
J d^/c 6fe exp(«A:a;), where V and he are the normalization 
volume and wavevector, respectively. The total energy 
density of the field at the present epoch, and in the 
absence of dissipation, is then given by spatial average, 
(...), which yields p% = {b^)/{STT) = l/lsnkl) J d^k\bk\^. 
Most primordial magnetic field models generate fields 
that are statistically isotropic and can be described by 
a power law spectrum, = BQ{k/kc)'^{n + 3)/47r 

for k < kc and zero otherwise, normalized such that 
(fo^) = Bq and Ic — 271 /kc corresponds to the field co- 
herence length. 

It is well known that injection of energy into the CMB 
between redhifts 10^ — 10^ and recombination, z ^ 10^, 
may yield observable distortions in the CMB blackbody 
spectrum, generically characterized by a chemical poten- 
tial ^ and a Compton y parameter jl^ . Given current 
upper limits on these parameters, i.e. < 9 x 10^^ 
and y < 1.5 x 10~^, one may infer that no more than 
~ 6 X 10~^ of the CMB energy density may have been 
generated in this redshift regime by nonthermal pro- 
cesses . Equilibration of photons to a blackbody dis- 
tribution has to proceed via photon-electron scattering, 
i.e. Compton scattering and double-Compton scatter- 
ing, slow processes due to the smallness of the baryon-to- 
photon ratio. Here, we are mainly concerned with black- 
body distortions resulting from diffusion (or mixing) of 
photons from different regions with different peculiar ve- 
locities. This process initially yields a spectrum which 
may be described by a superposition of blackbodies of 
different temperature, i.e. a Compton y distortions [p^ . 
When the mixing occurs at high redshift 2 ^ 4 x 10"' the 
optical depth of photons towards Compton scattering is 
still large, redistribution of photons in frequency space 
is possible, significantly diminishing the amplitude of the 
final y-distortion. Nevertheless, for redshifts z ^ 2 x 10^ 
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the final distribution will not contain the same photon 
number as a blackbody at the same temperature T, due 
to the inefficiency of photon number creating processes 
at large frequencies, of which double-Compton scatter- 
ing is the most rapid one. The resultant distribution will 
then be of the Bose-Einstein type at large frequencies 
with non-vanishing chemical potential /i. In contrast, 
when mixing occurs at lower redshifts z ^ A x 10^ spec- 
tral distortions are simply described by a non-vanishing 
Compton y parameter. For a more detailed discussion 
on the CMB equilibration process we refer the reader to 
Refs. jl8| and p9j, which have discussed generation of /i 



and y distortions resulting from primordial turbulence, 
very similar to the processes under consideration here. 

It has been recently suggested [|o| that the presence of 
pre-recombination magnetic fields affects the emission of 
gyroradiation by electrons into the CMB. However, the 
authors erroneously argued that this emission leads to 
large chemical potential distortions, even for very small 
magnetic fields, due to the efficiency of creating "extra" 
photon number at the fundamental, and small, frequency 
of gyration, hvg/T ^ 10-^(Bo/10-®Gauss)(z/10*^). Un- 
fortunately, the authors did not consider the inverse pro- 
cesses, i.e., the absorption of gyroradiation and inverse 
bremstrahlung, which even for very low redshifts are 
more rapid than the Hubble expansion due to the low 
frequency they occur at In spite of the emission 

of gyroradiation, the spectrum at low frequency remains 
Planckian due to the inverse reactions. 

We may now proceed to calculate the chemical poten- 
tial /i distortions resulting from magnetic field dissipation 
at high redshift. The evolution of chemical potential dis- 
tortions at large frequencies may be well approximated 

by ra 
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where Qb — dps/dt is the instantaneous energy dissipa- 
tion rate of the magnetic field. In this expression toe 
2.06 X lO^^s (1 - Yp/2)-^{nbh^)-^z' 
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characteristic time scale for double-Compton scattering, 
where Yp is the primordial helium mass fraction, f2(, is 
the fractional contribution of baryons to the critical den- 
sity, and h is the Hubble constant in units of 100 km s^^ 
Mpc^^. The solution to Eq. || at low redshifts is given 
by 

mW = 1.4/ dt^exp- - , (3) 

where — (4fo/5t£)c)^^^^ denotes the characteristic 
redshift for "freeze-out" from double-Compton scatter- 
ing. The quantity = 2.39 x 10^^ s determines the time- 
redshift relationship via t = t^jz^ for cosmic times well 
before matter-radiation equality. The characteristic red- 
shift equals = 2.49 x 10^ for 9.^,h^ = 0.0125, and 
y„ = 0.24. 



The energy dissipation rate may be computed as the 
sum of the contribution from all different magnetic field 
Fourier modes 
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dt 

d^fc|6fe|2 (2Imw)exp( -2 / Inicjdt 
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Note that, by employing &, the decrease in magnetic 
field energy density due to the expansion of the universe, 
which is not of dissipative character, does not contribute 
to Eq. IJ. Similarly, the real part of the frequency in 
Eq. does not change the total energy density in mag- 
netic fields and fluid flow. When and Im lo are in- 
serted in Eq. one obtains 
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It remains to substitute Eq. |^ in Eq. ^ and perform the 
double-integral. An analytic solution may be obtained 
for fee S> k^ijzy^ with fc^) — (15?^e'^T/i/?o)^^^7 vahd when 
the cutoff in the spectrum is at much larger wavevec- 
tors than the characteristic wavevector damped by pho- 
ton viscosity at redshift (see below). The solution is 
then simply 



n2 / lO \ ("+3) 

^0 ( 3/2 
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with X = 1.4r(n/2 + 5/2)r(3n/5 -f 9/5)2-("+^)/2 
(6/5) (ri + 3) a numerical constant of order unity. For ex- 
ample, K = 2.09, 1.10, 0.77, 0.78 for magnetic field spec- 
tral indices n — 1,0,-1,-2, respectively. The scale 
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i)zf/'^ has a simple interpretation. It is the scale that is 
damped by one e-fold at redshift z^. For the above values 
of rihh^ and Yp the corresponding comoving wavelength 

is Ad = {2Tr) /{k^zft^^) = 395pc. Note that the ratio of 
magnetic field energy density and CMB energy density 
is given by B^/{8Trp°) = 10-^^(50/3.24 x lO-SGauss)^. 

The limit may be put into context by comparing 
the prediction of Eq. |^ with the present upper limit 
from COBE/FIRAS, |^| < 9 x 10"^ at (95% confidence 
level) [0 for a given realization of the magnetic field 
spectrum. For example, when the spectral cutoff is taken 
at /cc = 7r/rg°° with rg™ « 0.9pc(TQCD/100MeV)-i 
the approximate comoving QCD-horizon, as well as un- 
der the assumption that the magnetic field energy den- 
sity is in approximate equipartition with the radiation 
energy density (i.e. Bq = 3.24 x 10"^ Gauss), one de- 
rives — 1.35. Note that there are no constraints on 
such magnetic field spectra from BBN, since a significant 
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fraction of the initial magnetic field energy density is al- 
ready dissipated by neutrino viscosity before neutrino de- 
coupling pl| . For comparison, a scenario which produces 
uncorrelated magnetic dipole fields in different QCD hori- 
zon domains, has a super- horizon power spectrum tail 
with n = [Q, and is unconstrained. Nevertheless, mag- 
netic field generation mechanism with resultant spectral 
indices of n = — 1 Q and n = — 2 Q have also been pro- 
posed. Furthermore, it is possible that processes, such 
as inverse cascades in a turbulent primordial fluid, may 
significantly increase an initially small coherence length 
of magnetic fields ||^ . 

The derived limit may be restated in an alternative 
form, independent of the initial statistical magnetic field 
properties given by the unknown Ic and n. Primordial 
magnetic fields with comoving coherence length « 400pc 
cannot have fields of strength 3 x 10"* Gauss since they 
would yield chemical potential distortions /i ^ 10""'. This 
is due to a nonthermal energy injection Aps/p-y ^ 10"^ 
induced by magnetic energy dissipation at redshift z w 
2 X 10^. One may derive a similar limit for magnetic 
field dissipation at redshifts z ^ 4 x lO** which would yield 
Compton y distortions in the CMB. For a magnetic field 
strength of w 3 x 10"* Gauss, magnetic fields damp in 
the photon diffusion limit above redshift z ^ 2 x 10^ with 
maximum damping scale « 0.3Mpc 0. Some addi- 
tional damping of even slightly larger scales « 0.6Mpc 
may occur shortly before recombination in the radiation 
free-streaming limit. Since a magnetic field strength of 
So ~ 3 X 10"* Gauss corresponds to a fractional CMB 
energy density of w lO"'', magnetic fields on such coher- 
ence scales would produce y distortions above the obser- 
vational upper limit. 

In summary, we have shown that spatially varying 
magnetic fields existing prior to the epoch of recombi- 
nation, may lead to observable distortions in the CMB 
blackbody. Distortions may result from the injection 
of non-thermal energy into the CMB due to the dissi- 
pation of magnetic field induced peculiar flows in the 
baryon-photon fluid. Our main result is, that magnetic 
fields of strength ^3 x 10"* Gauss with comoving co- 
herence length « 400pc, may not be present at red- 
shifts z « 2 X 10^, since they would yield chemical po- 
tential distortions /^^lO"^. Similarly, magnetic fields 
of ^ 3 X 10"* Gauss on scales ^ 0.6 Mpc are not al- 
lowed since they would induce y distortions shortly be- 
fore recombination. Limits on the spectral index of a 
stochastic magnetic field distribution, depending on the 
total energy density in the field and its minimum coher- 
ence length, are also given. These limits arc complimen- 
tary to limits on homogeneous primordial magnetic field 
(^0 ^ 3x 10~^ Gauss) and limits on larger scales discussed 
above. Our limits are stronger than the present Faraday 
rotation upper limits on these scales which for ^ 400 pc 
is only Bq < 7 x 10"^ Gauss pi. 
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